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Abstract. A linear irreversible thermodynamic framework of heat conduc- 
tion in rigid conductors is introduced. The deviation from local equilibrium 
is characterized by a single internal variable and a current multiplier. A 
general constitutive evolution equation of the current density of the internal 
energy is derived by introducing linear relationship between the thermody- 
namic forces and fluxes. The well-known Fourier, Maxwell-Cattaneo-Vernotte, 
Guyer-Krumhansl, Jeffreys type and Green-Naghdi type equations of heat con- 
duction are obtained as special cases. The universal character of the approach 
is demonstrated by two examples. Solutions illustrating the properties of the 
equation with jump boundary conditions are given. 



1. Introduction 

The increasing importance of micro- and nanotechnology initiated an intensive 
research in heat conduction [TJ [5] . Experiments show deviations from the classical 
Fourier theory [3J Uj and there are several theoretical developments to understand 
the nature of the deviations H H EH EQ . 

The starting point of all these approaches is the balance of internal energy 

(1) pe + d i q i = 0, 

where p is the density, e is the specific internal energy, and q l is the conductive 
part of the current density of internal energy, i.e., the heat flow. We are dealing 
with rigid heat conductors, the dot denotes the partial time derivative, d l is the 
space derivative of the corresponding physical quantity, and index notation with the 
Einstein summation convention is applied. In the phenomenological generalizations 
of Fourier's law of classical heat conduction, the constitutive equation ^ is modified 
by additional terms. The most important modifications are the following: 
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Here, |2| is the classical Fourier's law [T^], ^ is the Maxwell-Cattaneo-Vernotte 
(MCV) equation E3EHE], Eq. @ is the Guyer-Krumhansl (GK) equation [TB], 
([5]) is known as the Jeffreys type or lagging heat equation [17] , and ^ leads to the 
Green-Naghdi (GN) equation [18]. The heat conduction coefficient is denoted by 
A, the relaxation time by r, and a%, a 2 and b 2 are further material parameters. 
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The origin — the motivation and derivation — of these equations is manifold. Ki- 
netic theory offers several different approaches for Fourier's law and for MCV equa- 
tion. For instance, moment series expansion of the Boltzmann equation results in 
Fourier's law in first order and in the MCV equation in second order jTHl [2D] ■ Also, 
the GK equation was originally derived with the help of special collision terms in 
the Boltzmann equation characterizing phonon-lattice interaction |16[ 12 1] . 

The phenomcnological theories are also diverse regarding the origin of the above 
equations [T71 [22J [7]. Fourier's law is the consequence of non-negative entropy 
production in classical irreversible thermodynamics [23], and the MCV equation 
is the result of the deviation from local equilibrium, characterized by an internal 
variable [24] • The heat flow is the well-established candidate of that internal variable 
to obtain compatibility with kinetic theory [25J [THl [2U] ■ 

The origin of weakly nonlocal extensions is a question in a phenomenological ap- 
proach, too. There are indications that the Guyer-Krumhansl equation is connected 
to nonclassical entropy current density [2j)l [27J [25J . The lagging heat equation ([5| 
was suggested according to an analogy to rheological models and its relation to 
thermodynamics is not clear [7J. One may obtain it by diverse simple mechanisms 
of temperature equilibration, like heat flow splitting ! 29] , or two-step relaxation 
[501 I3T] . The Jeffreys type heat conduction is free from several problematic aspects 
of the MCV equation [32]. At last, the theory of Green and Naghdi is peculiar. 
They have introduced a special scalar internal variable, the time derivative of which 
is the temperature, and have analysed the thermodynamic consequences of the de- 
viations from local equilibrium in a nonstandard way |18j . They have obtained a 
model of heat conduction that cannot be reduced easily to Fourier's law, and which 
may result in zero entropy production: heat conduction without dissipation. 

We can see from the previous concise and not complete survey that the derivation 
and the motivation of the weakly nonlocal generalizations of Fourier's theory are 
manifold. Most of them are considered valid models of physical phenomena if their 
microscopic derivation provides a clear mechanism of the modifications. In some 
important cases — like the Jeffreys type or lagging heat equation — the role of the 
second law is not clarified. All kinds of thermodynamically consistent justifications 
introduce a deviation from the local equilibrium entropy density. In case of nonlocal 
extensions, a deviation from the classical form the entropy current density is also 
considered. 

In this paper, we derive a unified model of the previous constitutive relations 
with the help of irreversible thermodynamics. We need only two simple and general 
assumptions: 

- there is a deviation from the equilibrium state, conveniently and generally 
characterized by a vectorial internal variable [33J [3H [35] , and 

- there is a deviation from the classical form of the entropy current, conve- 
niently and generally characterized by an arbitrary function that we will 
call current multiplier [27j . 

In a linear approximation, both the internal variable and the current multiplier can 
be completely eliminated and one obtains a general constitutive equation for the 
heat conduction. We will see that the restrictions from the second law of thermo- 
dynamics are nontrivial and reduce the number of independent coefficients. In 
this paper, we apply the simple heuristic constitutive theory of irreversible thermo- 
dynamics in order to consider the restrictions of the second law. A detailed second 
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law analysis of rigid heat conductors with the help of Liu procedure in [27 shows 
that current multipliers require second order weak nonlocality. 

There is an important advantage of a phenomenological thermodynamic ap- 
proach based only on general assumptions. As long as the general conditions re- 
garding the deviation from the local equilibrium are fulfilled by any microscopic or 
mesoscopic assumption regarding the mechanism and the structure of the devia- 
tion, the consequences will be the same: this is the universality of noncquilibrium 
thermodynamics. To demonstrate this property we show, that if the deviation 
from Fourier's law is characterised by a gradient of a scalar field, then the model 
equations reduce to the parabolic two-step model, and that, if the deviation from 
Fourier's law is characterised by a general vector field, then the vector field should 
fulfil an MCV equation in order to satisfy the thermodynamic constraints. 

2. Irreversible thermodynamics of heat conduction 

2.1. The entropy production. For modeling phenomena beyond local equilib- 
rium, we introduce a single vectorial internal variable denoted by £*, and through 
this paper we consider rigid isotropic materials. The second law is given in the 
following form: 

(7) ps + d 1 ! 1 = a > 0, 

where s is the specific entropy, J 1 is the conductive current density of entropy and 
a is the entropy production. Then we introduce two basic constitutive hypotheses, 
which will be crucial in the following: 

(1) We assume that nonequilibrium entropy depends on the internal variable 
£ l quadratically: 

(8) »(e,a = «(e)-^ J . 

Here, m is a scalar material coefficient, sometimes called thermodynamic 
inductivity 25j. The quadratic dependence, with m = m(e,£ J ) can be 
justified by the Morse lemma and the requirement of entropy maximum at 
the nonequilibrium part of the basic state space, spanned by £ l [33J . m > 
because of the concavity of entropy. If m is constant then we obtain the 
following partial derivatives: 



(9) 



ds 



= -mC, 



where T is the equilibrium temperature. The thermodynamic relations are 
conveniently expressed by the following Gibbs relation: 

(10) de = Tds + mCTdC. 

In what follows, we make the assumption m = const., which is convenient 
in the following calculations, although it restricts the applicability of the 
theory. Then, the temperature is the function of the internal energy only 
and does not depend on the internal variable, as one can see from the first 
equation of 

(2) Our second assumption requires the generalization of the classical conduc- 
tive current density of the entropy. We demand that, without energy cur- 
rent, there is no entropy current in pure heat conduction. Therefore, in 
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the following we assume that, instead of the classical J 1 = q l /T form, the 
current density of the entropy is 

(11) ./' /ry. 

where B lJ is a constitutive function that will be constrained by the second 
law. This assumption was first introduced by Nyrri |36j . and applied for 
heat conduction in |26| . where the generalized multiplier of the current 
density of the internal energy was called current intensity factor. The recent 
terminology ('current multiplier') was suggested in |27j . 

The above two conditions lead to the following form of the entropy production: 

ps + d\P = -^d l q l - pmCC + d' l (B l3 q j ) = 

(12) (/)•'//" ) + <Vy ^B ij - ~S ij j - pmf f > 0. 

This formula is a quadratic expression of thermodynamic fluxes and forces, where 
the following characterisation is applied: 
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Table 1. Thermodynamic fluxes and forces 

The heuristic method of classical irreversible thermodynamics introduces a func- 
tional relationship between these quantities, assuming that the fluxes depend on 
the forces. Then in simple cases and for twice differentiable functions the mean 
value theorem of Lagrange ensures the validity of linear relations |37j . However, 
in our case the bilinear expression above does not permit such a simple character- 
ization, the natural choice of the constitutive state space spanned by the internal 
energy e and the vectorial internal variable £ l and their gradients is seemingly not 
without problems. However, in [27 j it was shown that, in case of a second order 
weakly nonlocal state space, the above choice of the entropy flux in the form of 
([2]) does not lead to a contradiction and, considering also that composite functions 
can appear in constitutive relations (see e.g. [35]), the consequent simple heuristic 
exploitation can be corroborated. 



2.2. Linear relations. Now we introduce such linear relationships that can be 
easily transformed into a convenient form. In isotropic continua the entropy pro- 
duction is an isotropic function of the forces, therefore, the most general linear 
relationship [3J5] between them introduces seven material parameters: 

(13) q l = h&BV-l 12 e, 

(14) mpC = l 21 d j B ij /.i". 

(15) B ij - ^ = k 1 d i q j + k 2 8 3 q l + k 3 d k q k S l: >. 

Here, l\, l\ 2 , I21, h, ^1, k 2 , k% are the isotropic scalar conductivity coefficients, and 
6 lJ is the Kronecker symbol. The non-negative entropy production requires the 



following inequalities for the material parameters: 

h > 0, h > 0, h > 0, k 2 > 0, 
1 



k 3 > 0, 



(16) 



L = hh 



:(h2 + l2ir>0. 



We do not assume Onsager or Casimir type reciprocal relations for the vectorial 
thermodynamic interactions characterized by the last two terms of the entropy 
production. Our internal variable may be a function of several of microscopic 
variables of different nature regarding time reversibility, some a type and others /3. 
This general approach was shown to be fruitful in formulating the thermodynamic 
framework of generalized continuum mechanics [401 141 j . 



We can eliminate the current multiplier from (13 1 and (14) with the help of 



(15 1. Moreover, after some simple steps, we can eliminate the internal variable 



from Eq.(13) and Eq.(14), too. If h ^ and m, l±, I12, &i> k 2 , k 3 are constants then 



we obtain the following constitutive relationship of the derivatives of temperature 
and of energy current density: 



(17) 



d , . 
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ot at 



Here, we have denoted the partial time derivative by d/dt and introduced the 
shorthands 

^12^21 



mp 



(18) 



Ai 

ai = Ai(fci + k 3 ), 
h = A 2 (fci + fc 3 ), 



h 



A 2 

a 2 = Aifc 2 , 
b 2 = X 2 k 2 . 



rap- 



We can see that (17 1 contains only five independent material parameter, as k\ and 
^3, l\ 2 and l 2 \ cannot be distinguished. Every coefficient in ( 18 ) is positive according 



to (|16j). The heat conduction coefficient deserves a special attention because 
(19) Wh - h 2 l 2 i = hh - l 2 s +l 2 a >0, 

where l s = (/i2 + ^2i)/2 and l a = (Z12 — fei)/2 are the symmetric and anti symmetric 



parts of the matrix in (13)-(14). As a consequence if Ai = 0, then A2 = follows. 



2.3. Special cases. We derive the following important special cases: 

(1) Fourier. If k\ — k 2 = k 3 = and I12 — 0, then directly from (13)-(15) we 
obtain Fourier's law, in the following form: 

(20) q i — Xidi-^ = —XdiT, 

where A = X1/T 2 = h/T 2 is the Fourier heat conduction coefficient. This 
is not apparent from ( 17 1 because ( 18 1 and ( 16 ) exclude the straightforward 
choice t — 0, A2 = 0, a\ — a 2 — 0, b\ — b 2 — 0. 

(2) Maxwell-Cattaneo-Vernotte. It is frequently mentioned 19J that extended 
irreversible thermodynamics, and the MCV equation, arises as a special 
choice of a vectorial internal variable as the conductive current density of 
the internal energy = q l [331135]. However, in our framework, the choice 
£ l = q l does not automatically lead to the MCV equation. In fact, (17) 
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shows that the choice of A2 = 0, a% = 0,2 = 0, b\ = hi = 0, leads to equation 
([3]). Thence ^i=0, and therefore, the internal variable is proportional to the 
heat flow according to (13). Ai = l^/h because of the last inequality of 
(fl6|). Since in (13) the first term of the right hand side is responsible for 



(21) 



(4) 



(5) 



(22) 



the classical Fourier type heat conduction, we may say he MCV equation is 
obtained if heat conduction is dominated by a Casimir type antisymmetric 
cross effect. 

(3) Jeffreys type. If a± = a% = 0, b\ = &2 = 0, we obtain the thermodynamic 
version of the Jeffreys type equation in the following form: 



X 1 d l - + X 2 ^- ( d l - 
T at \ T 



If Zi 7^ then A2 7^ follows, and the MCV equation is completed to a 
Jeffreys type equation. The emerging nonlinearity cannot be circumvented 
by assuming temperature dependent coefficients. 

Guyer-Krumhansl. If A2 = 0, b\ = 62 = 0, and Ai = AT 2 , then the GK 
equation Q is obtained. As l\ — is the simplest condition for A2 = 0, we 
see that the GK equation requires a Casimir type coupling of the terms in 



(13H14I 



General Green-Naghdi type. A GN type equation can be obtained if l 2 = 0. 
In this case, the Casimir type reciprocity l a = l\2 = —hi is the consequence 



of the last inequality of (16), and we obtain 



T m q 



Aid 1 



ot at 



Here, we have introduced similar notation of the coefficients as in ( 17), but 



they are combinations of the thermodynamic parameters different from that 



was given in ( 18 1 



(23) 



r = mp, Ai = Z Q , 
ai = Ai(fci + fc 3 ), 
h = A 2 (fci + fc 3 ), 



A 2 = mph, 
a 2 = \\k 2 , 
b 2 = X 2 k 2 . 



(6) Green-Naghdi type. A pure GN type equation requires 1% = 0, l\ = 0, and 
the Casimir type reciprocity l a = l\2 — Then we obtain 

(24) r ^ = Ai5*i + ai0*V + a 2 d"q\ 

The GN type equation may be nondissipative with zero entropy production 
when the entropy current density is classical, i.e. ki = It2 — k% = 0. 
An observation common to all the above cases is that heat conduction coefficient 
Ai is always nonnegative. 

3. Macroscopic universality 

In order to obtain an impression of the specific mechanisms that may lead to 
weakly nonlocal deviation from the classical entropy current, we consider here a 
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simple form of the GK equation, where, in addition to the conditions A2 = 0, 
b\ = b 2 = 0, we introduce k x = k 2 = 0. In this case, (17) simplifies to 



(25) rtf + q l = Xx&B, 
with B = l/T+ k 3 d k q k , t = mp/h and X x = l 2 a /l 2 . 

3.1. Two heat fluxes. As a first example, let us suppose that the deviation from 
Fourier's heat conduction is characterized by a vector field <f : 

(26) <f := q l + X s d% 

with Ag being a constant parameter. We may substitute q l into the simplified GK 



equation (25) and obtain the following condition for the unknown vector field: 



(27) r<f + <f = (A x - X s T 2 )d 1 ^ + (tX s - X 1 k 3 pc)d i f. 

Here we applied ([I]), together with the equation of state e — cT with constant 
specific heat c. Hence, with the choice of ^3 = rXs/{Xipc), our vector field <f fulfils 
an MCV equation with the A = A5 — A1/T 2 Fourier heat conduction coefficient. 
Therefore, the total heat flow q % is divided into two parts: q % — q l fulfils Fourier's 
law and q % an MCV equation. Apparently, there are two different channels of heat 
conduction characterized by a Fourier's law and an MCV equation, respectively 

CD. 



3.2. Two temperatures. We can introduce a different characterization of the 
deviation from Fourier's law by a scalar field T 2 in the following way: 

(28) 0T 2 = q* + X T d l T, 

where At and j3 are constant coefficients. Substituting q l to the simplified GK 



equation ( 25 ) again, one obtains the following condition: 



(29) d l (p{rf 2 +T 2 - T)) = (Ax + (A T - /3)T 2 ) 0*1 + (r A T - X^pc^f. 

Therefore, with the choice of k 3 — tXt/ (Xipc), Ai = ifi — Xt)T 2 one can find 
the condition that our scalar field T 2 fulfils a heat exchange equation. Apparently, 
there are two different temperatures that equilibrate by heat exchange. This is a 
well-known mechanism of heat conduction in complex materials where one of the 
components can be thermally excited independently of the other, e.g. in case of 
metals with different electron and lattice temperatures [30j [3Tj . 

Therefore we can see that both heat splitting and heat conduction of two- 
temperature materials are special cases of weakly nonlocal heat conduction and 
our approach with a general internal variable can be well interpreted in terms 



of the general weakly nonlocal heat conduction law (17). The current multiplier 
= BS l J is connected to the coefficients of the specific relaxation mechanism. 
These two examples show the real advantage of our weakly nonlocal generalization 
of heat conduction: it is based on minimal and general assumptions regarding the 
current density of the entropy and therefore the results are robust and universal. 

A straightforward microscopic interpretation of the internal variable could be a 
moment of the phonon distribution function. That view is suggested by the qua- 
dratic dependence of the entropy density of £. However, a detailed kinetic interpre- 
tation of current multipliers in the generalized entropy flux is an open question in 
the light of the above mentioned macroscopic ones. 

7 



4. Solutions of the generalized heat conduction equation 



In this section we solve a special form of the above equations for a one dimensional 
space variable x. Assuming a constant specific heat, e = cT we can eliminate the 



heat flux from Eq.(17l with the help of the energy balance and obtain 



\ m > + "57 ~ ~ — 7772 ^ ~ EZ^Im + *lK~ 2a . + X 2 k- 



dt 2 dt pcdx 2 T pcdx 2 dtT 1 dx 2 dt z dx 2 dt 2 

where k = k\ + k 2 + fc 3 . We can see that, due to the appearance of the reciprocal 
temperature, the equation is nonlinear and the Jeffreys type and GK equations do 
not coincide. In this equation, we have introduced constant conduction coefficients 



based on ( 13 )— ( 15 ) as it is natural in a thermodynamic approach. We can observe 
that, in this case, the Fourier heat conduction coefficient is temperature dependent, 
Ap = X1/T 2 7^ const. . 

In the following, we show some characteristic solutions of the differential equation 



(30 1 in case of jump boundary conditions. Initially we assume that a rod with length 
L has a fixed temperature T(0, x) = Tj and also <9tT(0, x) = 0. The far end of the 
rod is kept at this initial temperature, therefore, T(t, L) =Ti. The other boundary 
has the same temperature at the beginning, but after a short time the temperature 
jumps to a higher value Tf. For the sake of calculation convenience, we approximate 
the jump with the following smooth function: T(t, 0) = T J + (T / -T,)(l-e~ t / t i"»p) 2 . 

The used parameters are those of the NaF crystal experiments of Jackson and 
Walker [33J H3I- Therefore the specific heat is c = 2.774W/Kkg, the density p = 
2866kg/m 3 , the Fourier heat conduction coefficient Xp = 20500W/mK, therefore, 
Ai = 4.61 x 10 6 WK/m, the length of the rod is L = 8.2pm, and the parameters of 
the initial condition are Tj = 15K, Tj — Tf = 0.1K, ij ump = O.bps. 

Figure [T] shows the solution of the Fourier equation. The linear stationary solu- 
tion T(x) = Tj + Tf x is reached at the end. 

In the Figure [2] a solution of the Maxwell-Cattaneo equation is shown whith 
t = 0.68/is and the other parameters are the same, A2 and k are zero. The relaxation 
time is extracted from the measured v c — 1.95mm//xs. Here one may observe 
a ridge due to the finite propagation speed, and the damping of the ripples and 
of the reflected pulse due to the dissipative effect of heat conduction. Figures [3j 
and [4j show the solutions of the Jeffreys type and the Guyer-Krumhansl equations, 
where r is the same as in Figure [2] and A 2 = 0.05JK/m (Jeffreys case) and k = 
10~ 12 ms 3 /kgK (Guyer-Krumhansl case) respectively. In both cases, the relaxation 
to the stationary solution is faster. In Figure [5] one can see a solution where none 



of the coefficients of equation (30) is zero . There the coefficients are r = 0.68/is, 
A2 = 0.05JK/m, k = 10~ 12 ms a /kgK. In this case, one may observe the change of 
the speed of the signal propagation. 

The solution of the nondissipative GN type equation is given in Figure [6] with 
the same coefficient values as in case of MCV equation, r = 0.68/is, Ai = 4.61 x 
10 6 WK/m. Figures [7] and [8] show that Jeffreys or Guyer-Krumhansl dissipation 
do not damp the reflected travelling wave (parameter values: r = 0.68/is, A2 = 
0.05JK/m, respectively k = 10~ 12 ms 3 /kgK in Figure [8j Finally Figure [9] shows 
that, for the general GN equation, the damping effects is due to the additional 
terms are different, from which can be observed for the GK- Jeffreys type one. 
Using the t = 0.68^s, A 2 = O.lJK/m, k = 2 x 10~ n ms 3 /kgK parameter values, 



Figure 1. Fourier equation, r = 0/j.s. 



15.1 



T [K 




x [ram] 



Figure 2. MCV equation, t = 6.8^s. 



one gets a qualitatively different damping than with the same parameter values for 



the heat conduction equation ( 22 ) 



5. Summary and conclusions 

We have introduced a vectorial internal variable and a current multiplier as a 
simple and general characterization of the deviation from local equilibrium of the 
entropy and the entropy current density functions, respectively. Then we have 
derived the resulting connection between the heat flow and temperature, and ob- 
tained a common generalization of several known heat conduction equations. Con- 



sequently, the emerging heat conduction equation (17) is independent of the micro- 



Figure 3. Jeffreys type equation, r = 6.8/xs, A2 = 0.05JK/m. 




Figure 4. Guyer-Krumhansl equation, r = 6.8/zs, k = 10 12 ms 3 /kgK . 



or mesoscopic mechanism of heat conduction as long as the general conditions of 
this approach — namely, the second law and the existence of extra field quantities — 
are fulfilled. In this sense, the found heat conduction equation is universal. We 
have demonstrated this kind of universality by the derivation of two common simple 
heat conduction mechanisms in our framework. 

The analysis of the different specific examples in the light of the thermodynamic 
constraints, stemming from the non-negative entropy production, revealed that: 

If the main Fourier coefficient (li) is not zero then MCV is always extended 
by the characteristic nonlocal term of the Jeffreys type equation. 
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Figure 5. Jeffreys-Guyer-Krumhansl equation, r = 6.8^s, A 2 = 
0.05Jif/m, k = 1Q- I2 ms 3 /kgK. 




Figure 6. Nondissipative Green-Naghdi type equation r = 



Pure MCV, GK and GN type equations are related to a Casimir type 
cross effect between the internal variable and thermal parts of the entropy 
production. 

The nondissipative wave equation of the GN type model was derived and 
interpreted in the framework of non-equilibrium thermodynamics. 
The strictly linear thermodynamic constitutive equations < 13 )— ( 15 1 with 
constant coefficients lead to a nonlinear evolution equation of the heat flow, 
where the nonlinearity cannot be absorbed in special temperature depen- 
dent coefficients. 



Figure 7. Green-Naghdi type equation with Jeffreys damping 
t = 6.8^s, A 2 = OmJK/m. 




The general heat conduction equation (17 1 is remarkably stable in numer- 
ical calculations, whenever the inequalities ( 16 ) are fulfilled. (For the link 
between stability and the second law of thermodynamic, see e.g. [IS]) 
With the natural choice of the internal variable as the non-Fourier part of 
the heat flow, a physical interpretation of our internal variable was given. 
We found that two-temperature relaxation and heat flow splitting models 
can be interpreted in terms of the current multiplier. 



T[K] 



15.1 



1^ 
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Figure 9. Green-Naghdi with both kinds of damping, r = 6.8/is, 
A2 = O.lJK/m, k = 2 x 10~ n ms 3 /kgK, longer duration. 

It is worth mentioning that ([2| trivially corresponds to the usual generalisation 
of the entropy current density, when an additive constitutive term, the so-called 
Miiller's K l vector is introduced (see e.g. [46], p441). 

Another comparison that one can perform is related to the GENERIC approach 
in [5]. Equations (42)- (49) in [5] are in principle more general than ours from several 
points of view. On the other hand, the structure and derivation of the equations 
is rather different, therefore, the comparison is not straightforward. For example, 
the suggested evolution equations are balance like and the particular form of the 
entropy flux (first part of Eq. (46) in [9) looks structurally different from our 
simple generalization, preventing a direct comparison. 
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